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Abstract 

We study a two loop induced radiative neutrino model with global ?7(1) symmetry at 1 GeV 
scale, in which we explain the X-ray line at 7.1 keV of dark matter candidate reported by XMN- 
Newton X-ray observatory using data of various galaxy clusters and Andromeda galaxy. We also 
discuss Higgs sector, lepton flavor violation processes, and a physical Goldstone boson. 

Keywords: Radiative seesaw, Goldstone Boson, X-ray Line 


"Electronic address: hokada@kias.re.kr 


1 




Lepton Fields 

Scalar Fields 


Ll 


Nr 

$ 

xt 

X2 

+ 

SU{2)l 

2 

1 

1 

2 

1 

1 

1 

U{1)y 

-1/2 

-1 

0 

+1/2 

+1 

+ 1 

0 

U{1) 

— X 

— X 

x/3 

0 

2x 

2x/3 

-2x/3 


TABLE I: Contents of lepton and scalar fields and their charge assignment under SU {2)lxU (l)y x 
C/(l). 

I. INTRODUCTION 


One of the promising scenarios simultaneously to explain between Neutrinos and dark 


matter (DM) physics is to generate neutrino masses at multi-loop level 
in which DM could be a mediated held in the neutrino loop. 

As for the DM sector, X-ray line signal at 3.55 keV from the analysis of XA 


47|, |48l-l63|, 


m 


X-ray observatory data of various galaxy clusters and Andromeda galaxy 




N-Newton 


66] can be 


understood by decaying scenario, in which the DM mass should be 7.1 keV and mixing 
angle between DM and the active neutrinos sho uld b e sin^ 29 ~ 10“^°. Due to these simple 
implications, many works have been studied 67Hll8l |. ^ 

In our paper, we propose a two loop induced radiative neutrino model with a global U{1) 
symmetry, in which such a small mixing between DM and neutrinos can be generated at 
one-loop level. Since the DM mass is expected to be very tiny, it does not link to the neutrino 
masses. Observed relic density can be thermally obtained by the annihilation process with 
a Goldstone boson (GB) pair that is a consequence of the global U{1) symmetry. 

This paper is organized as follows. In Sec. II, we show our model building including Higgs 
sector, neutrino masses, and lepton havor violations (LFVs). In Sec. Ill, we analyze the DM 
properties including relic density and X-ray line. In Sec. IV, we numerically analyze the 
allowed region to satisfy all the conditions. We conclude and discuss in Sec. V. 


^ On the other hand, there are refuting papers 119lll2l| that the 3.55 keV line may imply atomic transitions 
in helium-like potassium and chlorine are more likely to be the emitters. 
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II. MODEL SETUP 


We discuss a two-loop induced radiative neutrino model. The particle contents and 
their charges are shown in Tab. [H We add three gauge singlet Majorana fermions Nji, two 
singly-charged singlet scalars xt)y a neutral singlet scalar ip to the SM. We assume 
that only the SM-like Higgs $ and (p have vacuum expectation values (VEVs), which are 
symbolized by v/y/2 and v'/\/2 respectively. We also introduce a global 17(1) symmetry, 
under which a; 7 ^ 0 is an arbitrary number of the charge under the 17(1) symmetry, and 
the neutrino masses are induced at the two loop level after the spontaneous 17(1) symmetry 
breaking by 2x as shown in Figure [2] (as well as Figure [5]). Notice here that SM-like Higgs $ 
should be neutral under the 17(1) symmetry not to couple quarks to Goldstone boson through 
chiral anomaly to be consistent with the axion particle search. Otherwise its breaking scale 
should be very large( that is greater than 10^^ GeV) , which is inconsistent with the observed 
value V ~ 246 GeV. 

The renormalizable relevant Lagrangian and Higgs potential under these symmetries are 
given by 


-Cy - {yi)ij{LL)i^{eR)j 


■ LL^xt + 9ijNR,e%^X2 + ^{yN)iMN‘^)i{NR)j + h.c.. 


(H.l) 


V 





+ |2 
2 I 


+ -I- Ao((^^Xi'X 2 + h-c.) -I- -I- Xx2\xt\‘^ 


+A.iv^r+A<i>id>r++ a,,^ix^ n<i>p 

+ Ax.,,|x^l Vr + Ax.$|x^n$r + A^$|9:^ri<hr, (H.2) 


where i = l — 3, j = l — 3, the hrst term of Cy can generates the SM charged-lepton masses 
me = yev/V^ after the electroweak spontaneous breaking of $, and we assume Aq to be 
real. The Majorana mass (Mat = yNv'/ \/2) is generated after the spontaneous breaking of 
p>. The scalar helds can be parameterized as 


$ = 


w ' 

v+4>+iz 

72 


ip = 


v' + a 

TT' 


,iG/v' 


(11.3) 


where v ~ 246 GeV is VEV of the Higgs doublet, and and are respectively (non¬ 
physical) GB which are absorbed by the longitudinal component of W and Z boson. Since 
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the CP even bosons (0, a) and singly-charged bosons {xt^xt) each other through the 
term and (p‘^XiX 2 respectively, each of the resulting mass eigenstate and mixing 

matrix is reparameterized as 60 1 


a 

^ vt 

hi 

5 

xt 

= 0^ 

ht 
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h2 


xt 


ht 


(11.4) 


where h 2 is the SM-like Higgs, hi is an additional CP-even Higgs mass eigenstate, {hf, h^) is 
the singly charged boson mass eigenstate, and each of O and H is 2 x 2 unitary mixing matrix. 


The lower bound of t 


80 GeV < 


le singly-charged boson can be obtained by the LEP experiment; he.. 


A. Invisible decay of the SM-like Higgs (/12) 


The current experiment at LHC suggests that the invisible branchin g ratio of the SM 
Higgs (Hinv) is conservatively estimated to be less than 20 % jl23l . Il24| . There are three 
invisible modes: ^2 —?■ 2G, ^2 —?■ 2Nk, and ^2 —)■ 2hi, and their decay rates (Pinv) are given 
by 


Tinv = r(h2 —)■ 2G) P(h2 —)■ 2NjS) -\- —)■ 2hi), 

nir. ^ 2G) = t|lM, 

rrih^Mlr |Pi2p , 

P(h 2 ^ 2Nk) = ' ( 1 - 


P(h 2 —)■ 2hi) — 


167rn'^ 

\lJ'h2hih\ I 


AM\ 


2 \ 3/2 
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h2 


Ami 


hi 


2 ’ 


(11.5) 

( 11 . 6 ) 

(11.7) 

( 11 . 8 ) 


?>2'Kmh2 V ^h2 

hh2hihi = yii{^ip^y22V + 6A^V2inO + 2A<^$l/iiHi2(V2in -I- V 22 V') + 4“ 6A$V22'r'). 

(11.9) 

Then we have to satisfy the following relation 


P < 

-1- inv \ 


Anv|H2-'2 


I-B 


221 pSM 

Higgs) 


( 11 , 10 ) 


where Pp^gg ~ 4 MeV is the total decay width of the Higgs boson at = 125 GeV. 
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FIG. 1: One-loop contributions for the pure quartic couplings of charged bosons where 

a = 1,2 in the left-hand side. 


B. Constraint on the charged bosons 


The vacuum stability should be satished for the pure quartic couplings of charged bosons; 

model, there exist some loop contributions to these couplings, by mediating 
bosons and fermions. Here we consider this issue at the one-loop level. As a general state¬ 
ment, we obtain a negative contribution relative to the tree-level coupling in case of boson- 
loop. On the other hand, a relatively positive contribution is given in case of fermion-loop. 
We have only a boson-loop contribution to the A^^ coupling as shown in the left-hand side 
of Figure [U on the other hand, we have both contributions to the A ^2 coupling as shown in 
the right-hand side of Figured! Then each of the condition up to the one-loop level can be 
given as 


\ one—loop ^ \ _ 3 

^ 32(47r)2 (M2^ - M|)3 _ 


{Ml + Ml) In 


M\ 


2 1 

Xa 


Mi 


- - Ml) 


AMpf.Misf ^ \\ f dadbdcdd6{a + b + c + d — 1) 

2^ (47r)2 ^dkjdikgHg.j) J 0^2^ + 5^2 + + dMl 


i,j,k,l 


> 0 , 


( 11 . 11 ) 


^Xa = ^Xa + = 771$ + = -A$n^ (11.12) 


where a = 1, 2, we have used the tadpole condition of $ for the last term, and we neglect 
terms that are proportional to or {vv')"^ because v'{k, 0{1) GeV) is very smaller than 
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FIG. 2: Dominant contribution to the neutrino masses at the two-loop level. 


the mass scale of (that is 0(100) GeV).^ The second term of the right-hand side in 

Eg. (III.lip corresponds to the contribntion of the boson-loop diagram as shown in the left- 
hand side of Figure [H and The third term of the right-hand side in Eg. (111.lip corresponds 
to the contribution of the fermion-loop diagram as shown in the right-hand side of Figure [H 

C. Neutrino mass matrix 

At first we redefine relevant terms in terms of the mass eigenstate as 

C-y ~ ■ U,ht + 9.j{0‘')^Nr:Ah- + + h.c. 

= /“Li, ■ + h.c., (11.13) 

where /“. = /p(Ol)ia, = fi'p(Ol) 2 a, and yj^Vx^l^ = y\,. 

Then the dominant contribution to the active neutrino mass matrix mp is given at two- 
loop level as shown in Figure [21 and its formula is given by 

("ii.)p = (11.14) 

^ 3 

^ The perturbativity and avoiding the global minimum can be straightforwardly satisfied if each of quartic 
coupling does not exceeds Itt and each of the sum of mass terms and the couplings are positive. 
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where Caj = {m(,Jmy^+Y, Mn^. = Mn 2 , M^^), and me = {me,m^,mr). Notice here 

that another contribution to the neutrino mass matrix (Zee-Babu-like diagram as shown 
in Figure Ej) can be tiny enough to be neglected in our case. This formula is found in the 
Appendix. 

The obse rved mixing matrix; PMNS(Pontecorvo-Maki-Nakagawa-Sakata) matrix 
(h^PMNs) |l25l |. can always be realized by introducing the Casas-Ibarra parametrization |l26| . 
where it is given by 




1-3 1,2 
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(11.16) 


Faj = O is an arbitrary orthogonal matrix with complex values. Then the 

neutrino mass eigenvalues can be given by 


mlm^, = t/pMNS 


/ 2 

0 
0 


m 


V2 


\ 


0 m 


U. 


t 


PMNS- 


V3 


( 11 , 17 ) 


When m^+ = 0(100 GeV), F k, 0.01 is obtained. Then written in terms of 

masses of neutrinos, charged-leptons, and right-handed neutrinos as 





( 11 , 18 ) 


Depending on the charged-lepton masses, we can estimate a typical order of the Yukawa 
coupling as 


= 0(10-7-10-3), (11.19) 

a 

where we fix = 0.01 eV, and M^v = 0(10-® — 10-^) GeV. 
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D. 


Lepton Flavor Violations (LFVs) and lepton universalities 


process 


The constraints from the —)■ £t£f.£i process |l28| can be given at one-loop level as 


{dhSa)ji{9o.9b)ik T i9o.9b)jki9b9a)ii 


256 


{flfa)Aflfb)ik + {flfa)jk{flfb)i£ 


8(47r)2(m^± — m^±) 


In 


'^a 

ml± 


< 


a 


ijki 


( 11 . 20 ) 


[GeV]2’ 

where each of Cijki is given by C'^.eee ~ 2.3 x 10"^ Creee ~ 0.009, Cree^L ~ 0.005, 


0.007, Cruee ~ 0.007, C. 


Tfiefi 


^ 0.007, and 0.008 


1291]. Notice hereafter that we 


In 


assume to be << m^± . Notice here that the form —reduces to — in 

... m.^ I —m. I m. , 


1 ( 2 ) 


m ,—m- _i_ 
hn 


the limit a{b) —)■ b{a). 


‘6(a) 


2. i, 7 process 


The constraints from the £ • 7 process |l30| can be given at one-loop level as 


{9i9a) 

Arn^. 


ml^ 

% 


< 


[GeV]4’ 


where each of Cij is given by G^e ~ 1-6 x 10 Cre ~ 0.52, Cr^ ~ 0.7 129 |. 


( 11 . 21 ) 


3. li /£j universality 


The constraint of the £il£j universality 13l| is given as 


n 


Il'n. 




jk 


m.± 


< 


a 


i/j 


[GeV]2’ 


( 11 . 22 ) 


wher e i ^ j 7 ^ fc, and each of Ct/j is given by G^/e ~ 0.024, Cr/^ ~ 0.035, and Cr/e 


0.04 


129|. 
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FIG. 3: The dominant contributions of the X —)■ i/lj process to explain the X-ray line. 


III. DARK MATTER 


We consider a fermionic DM candidate X[= iVi), which is assumed to be the lightest 
particle of iVj. And we focus on the explanation of the X-ray line at 3.55 keV, since X 
decays into active neutrinos and photon at the one-loop level after the spontaneous U{1) 
symmetry breaking as shown in Fgiure [3l Then the mass of DM Mx{= is hxed to 

be around 7.1 keV wit h a s mall value of the decay rate divided by Mx] he., 4.8 x 10“^® ;< 
~ ^ 10“'^® jin7 |. ® In our case, this form can be rewritten as 

1 O 1 Q . . 

jl f 1 3 — dCaj -|- 2 InCa^, 


4.8 X 10 -"® < 

~ 4(47r)4 


{fa)jk{9l) 

a j,k 


ml^ 

'I'a 


(1 


Ineaj)® 


< 4.6 X 10 -"®, 


(III.l) 


where CTem ~ 1/137 is the hne structure constant. 

Relic density. The observed relic density can be thermally generated through the process 
of the GB hnal state due to the global 77(1) symmetry. The relevant Lagrangian is given by 

Cdm = ^ V VrPhd^Gd^G). (III.2) 

6n V ^^ 

h=l,2 


Notice here that the term oi yx cannot contribute to the annihilation process of 2G hnal 
state, since this term does not connect to the other terms with GB interactions due to the 


absence of the Lorentz index of y. Then the relativistic cross section of X is given by 


2567rn'^(m^^ — s)^ 


(III.3) 


^ This bound is derived from sin^ 29 = {2 — 20) x 10 


104 |. 
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FIG. 4: Allowed scanning points in terms of Mn and r(X —)• v^)/Mx, where these points satisfy 
all the constraints of LFV processes and lepton universalities as discussed in Section IIIDl 


where we neglect t he c ontribution of the SM Higgs h 2 .Then the thermal averaged cross 


section is given by 


182 1 


(<7Wrel) 


ds^Js - Mj.aVre\{s)Ki{^/sx/Mx)x 


(III.4) 


16M|.[A:2(a:)]2 

where x = Mx/T at the temperature T, and each of it' 1^2 is the modihed Bessel function 
of the second kind of order 1 and 2. Finally the formula of the relic density [VLxh"^) is 
approximately given by 

1.07 X 10^ 


Vtxh^ 


•ya.(Xf)Mp, /“ [GeV] ’ 


(III,5) 


where the lower index of / represents the time of freeze-out, and each of Mpi 1.22 x 10^® 
GeV and g*{xf) is the Planck mass and the nu mber of the relativistic degrees of freedom. 
To obtain the observed relic density Vtxh'^ ~ 0.12 133| . we hnd a solution at n' = 1 GeV and 
0.051 GeV, where we use hxed reasonable values |Vii| = 1, Xf = 20, and g*{xf) = 100 
for brevity. We will adopt this solution in our numerical analysis. 
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IV. NUMERICAL ANALYSIS 


In this section, we search the allowed region for Yukawa couplings / and g. First of all, 
we £x the following relevant parameters for simplicity 


w' 1 GeV, ruh-^ ^ 0.0506 GeV, mh 2 ~ 125 GeV, m^± 152 GeV, 193 GeV, 

124 GeV, ^ 121 GeV, 


Mx-7.1 X 10"® GeV, V 


-1.00 1.54 X 10 


-3 


, O 


1.54 X 10-3 1.00 

A<i, ^ 0.130, 0.0984, ^ 0.254, Aq ~ 9.281, A^^ 0.962, A^2 ~ 0.614, 

A^^$ 0.864, A^2'1> ~ 0.745, A^^^ 0.0984, A^^^p ~ 0.915, 


-1.00 3.31 X 10-3 

3.31 X 10-3 1.00 


(IV.l) 


where the condition of vacuum stability in Eq. fill. lip is satisfied. Now we can check whether 
the above input parameters satisfy the constraint of the invisible decay of SM-like Higgs in 
Eq. flll.lOp . Our invisible decay width is computed as Finv 0.0464 that corresponds to the 
left-hand side of Eq. flll.lOp . On the other hand, the value of the right-hand side of Eq. flll.lOp 
is equal to 0.2. Hence these input parameters satisfy this constraint. 

Then the order estimation of Yukawa couplings in Eq. dll.lOp reduces to fg^ = 0(10-“^ — 
1), which can reproduce the sizable scale of neutrinos (0.01 eV). Within this estimation, we 
randomly select hfteen values of Yukawa couplings as 


Mx = Mx 2 (= Mx,) G [10-3,10-"] [GeV], (IV.2) 

/i,-e [5 X 10-3,0.1], (IV.3) 

r?,, e [5x10-3, 0.1] for ((1,1), (1,2), (1,3)), (IV.4) 

g'^g,, G [10-^10-^ for (^,j) = ((1,1), (1,2), (1,3)), (IV.5) 


and we take 1000 scanning sample points over these above ranges, where we expect that a 
tiny value of g' plays a role in explaining the X-ray line in order to satisfy Eq. fllll.ip . Figure 0] 
shows allowed scanning points(506) in terms of Mx and F(X —)■ i"'y)/Mx, where these points 
satisfy all the constraints of LFV processes and lepton universalities as discussed in Section 

KDl 
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V. CONCLUSIONS AND DISCUSSIONS 


We have constructed a two-loop induced neutrino model with a global U{1) symmetry, 
and we have shown allowed scanning points(506) in terms of and r(X u'y)/Mx in 
Figure m where these points satisfy all the constraints of the scale of neutrino masses, LFV 
processes, lepton universalities, relic density of DM (at v' = 1 GeV and nih^ 0.051 GeV), 
and X-ray line at 7.1 keV of the DM mass reported by XMN-Newton X-ray observatory 
using data of various galaxy clusters and Andromeda galaxy. Here DM does not link to the 
neutrino masses, since the DM mass is very tiny as well as g' coupling. 

As another phenomenological point of view, GB arises a discrepancy of the effective 
number of neutrino spec ies in the early Universe, which is denoted by AA'^ff. The recent 
data reported by Planck 133 | shows AA^"efr < O.OdlgJg at the 95 % conhdential level. In our 


case, AAgff is about 0.045, where we suppose that GB typically decouples from the plasma at 
temperatures just above the QGD phase transition(150-200 MeV) and our effective number 
of relativistic degrees of freedom is about 70. Therefore, our model can evade this constraint. 
It is also worth mentioning th at co smological issues such as an effect on cosmic microwave 
background via c osmi c string |l34| and stellar energy loss via the photoproduction process 
7 -|- e —)■ e-f GB |l35| may put some constraints only when it couples to an axial vector 
current. However GB in our model has a vector current only, therefore we can evade these 
constraints, too. 


Appendix 

The formula of the Zee-Babu type neutrino mass can be given by 


{jnv)nm — 

a,b j,k,l 


X / dxdydz5{x + y + z-l) [ dx'dy'dz' 


6{x' + y' + z' — 1) 


(l/2 - y)y'Xf^+ - {yXx^ + zX,^+)z‘ 


7, (V-1) 


where Xt = {nif/MY, and M = Max(m.+ , m.+ , Mat^) Max(m.+ , m.+), as shown in 
Figure [5l In our scanning range of the numerical analysis, the typical scale of the contribution 
to the neutrino masses is around 0(10“®) eV that is negligible. 
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FIG. 5: Subdominant contribution to the neutrino masses with Zee-Babu type. 
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